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Abstract Unique among neutron stars, IE 1207.4—5209 
is an X-ray pulsar with a spin period of 424 ms that con- 
tains at least two strong absorption features in its energy 
spectrum. This neutron star is positionally coincident 
with the supernova remnant PKS 1209—51/52 and has 
been identified as a member of the growing class of radio- 
quiet compact central objects in supernova remnants. 
From previous observations with Chandra and XMM- 
Newton, it has been found that the IE 1207.4-5209 
is not spinning down monotonically as is common for 
young, isolated pulsars. The spin frequency history re- 
quires either strong, frequent glitches, the presence of a 
fall-back disk, or a binary companion. Here, we report 
on a sequence of seven XMM-Newton observations of 
IE 1207.4—5209 performed during a 40 day window be- 
tween 2005 June 22 and July 31. Due to unanticipated 
variance in the phase measurements during the observa- 
tion period that was beyond the statistical uncertainties, 
we could not identify a unique phase-coherent timing so- 
lution. The three most probable timing solutions give fre- 
quency time derivatives of +0.9, —2.6, and +1.6 x 10~ 12 
Hz s _1 (listed in descending order of significance). We 
conclude that the local frequency derivative during our 
XMM-Newton observing campaign differs from the long- 
term spin-down rate by more than an order of magni- 
tude. This measurement effectively rules out glitch mod- 
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els for IE 1207.4—5209. If the long-term spin frequency 
variations are caused by timing noise, the strength of the 
timing noise in IE 1207.4—5209 is much stronger than 
in other pulsars with similar period derivatives. There- 
fore, it is highly unlikely that the spin variations are 
caused by the same physical process that causes timing 
noise in other isolated pulsars. The most plausible sce- 
nario for the observed spin irregularities is the presence 
of a binary companion to IE 1207.4—5209. We identi- 
fied a family of orbital solutions that are consistent with 
our phase-connected timing solution, archival frequency 
measurements, and constraints on the companions mass 
imposed by deep IR and optical observations. 
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1 Introduction 

There exist a handful of enigmatic X-ray point sources, 
very likely young neutron stars, positionally coincident 
with supernova remnants (SNRs) whose nature remains 
uncertain. These objects are commonly referred to as 
central compact objects (CCOs) that are characterized 
by soft, thermal X-ray spectra and an absence of ordi- 
nary pulsar activity such as radio pulsations, 7-ray emis- 
sion and pulsar wind nebulae (Pavlov et al. 2002a, 2004). 

The CCO IE 1207.4-5209 (1E1207 hereafter) in the 
PKS 1209—51/52 SNR is a particularly interesting mem- 
ber of this class in that it is both an X-ray pulsar (with 
a 0.424 s period; Zavlin et al. 2000) and the only CCO 
found to possess prominent absorption lines in its spec- 
trum (Sanwal et al. 2002). The X-ray energy spectrum is 
best modeled with a continuum blackbody component of 
temperature kT ss 0.14 keV and at least two broad ab- 
sorption lines centered at 0.7 and 1.4 keV. The strength 
of these lines depend upon the rotational phase of the 
pulsar (Mereghetti et al. 2002). Two additional features, 
at 2.1 and 2.8 keV, have been reported (Bignami et al. 
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2003); however, their validity has been questioned (Mori, 
Chonco & Hailey 2005). 

The physical origin for the spectral lines remains un- 
known. Sanwal et al. (2002) have concluded that these 
lines cannot be associated with transitions in Hydrogen 
atoms and argued that neither electron nor proton cy- 
clotron resonance could cause these features. These au- 
thors suggest that the lines could be due to absorption 
by once-ionized Helium in a magnetic field B ~ 2 x 10 14 
G (see also Pavlov & Bezchastnov 2005), while Hailey 
& Mori (2002) and Mori & Hailey (2006) argue that the 
lines could be formed in an Oxygen atmosphere with 
B ~ 10 n -10 12 G. As different interpretation imply very 
different magnetic field strengths, measuring the field 
strength of 1E1207 would be most important for under- 
standing the nature of the spectral lines. For isolated 
pulsars, the most straightforward method for estimating 
dipole magnetic field strengths is to measure the spin 
frequency v and its time derivative (spin-down rate) v. 

Zavlin, Pavlov & Sanwal (2004) studied the spin evo- 
lution of 1E1207 using a compilation of Chandra and 
XMM-Newton observations covering 3.5 years. They found 
that the spin frequency of this pulsar is not steadily de- 
creasing as one would expect for a magnetically-braking 
dipole. Instead, the frequency evolution was quite er- 
ratic, leading Zavlin et al. to consider three possible ex- 
planations: («) the star is undergoing frequent glitches, 
(ii) the star is surrounded by a debris disk that influ- 
ences its spin evolution through accretion and propeller 
torques, or (iii) the star is a member of a (non-accreting) 
binary system. Another possibility is that the spin-down 
of the star is influenced by timing noise, a ubiquitous 
property of isolated neutron stars whose physical origin 
is unclear. The sparse frequency history of 1E1207 could 
not distinguish between these models. 

Here, we report on a sequence of seven XMM-Newton 
observations of 1E1207 that were designed to perform 
phase-coherent timing in order to precisely measure the 
pulse frequency and frequency derivative of the source. 
A measurement of the local frequency derivative would 
help us to distinguish between the possible scenarios. 
Below, we describe the observation (§2) and our timing 
analysis of the XMM-Newton data set (§3), and discuss 
the resulting constraints on the physical mechanisms for 
the spin-frequency evolution in 1E1207 (§4). 



2 XMM-Newton observations 

During a 40 day interval between 2005 June 22 and July 
31, XMM-Newton observed 1E1207 seven times, with the 
EPIC PN camera as the primary instrument. The first 
three and last two pointings had effective exposures of 
10—15 ks, while the central (fourth) exposure was about 
45 ks. The fifth pointing of about 7 ks exposure was 
shorter than planned because of strong background con- 
tamination. The spacing between consecutive observa- 



tions was approximately 15, 5, 1, 1, 5, and 15 days. The 
spacing and durations of the XMM-Newton pointings 
were planned in such a way as to allow phase-coherent 
timing of the pulsar over the full time span of 40 days (see 
§3). The exact exposures, observing epochs, and other 
observational details for these observations are listed in 
Tabled 

For each observation, the PN camera was operated in 
small window mode, with 5.6 ms time resolution. Start- 
ing from the observation data files, all data were pro- 
cessed using XMMSAS version 6.5.0. After running the 
tool epchain, we extracted light curves from the ob- 
served field of view minus a circular region that included 
1E1207. These light curves were used to identify and 
filter out periods of high background. Source (plus back- 
ground) counts for timing analysis were extracted from 
a circular region for each observation and filtered using 
standard criteria and the good time intervals we deter- 
mined. The radii of the extraction regions, 35" or 20", are 
listed in Tabled A smaller radius of 20" was used to im- 
prove the signal-to-noise ratio for the three observations 
where the background rate within the good time intervals 
was elevated. The filtered event lists were barycentered 
to the location R. A. = 12 h 10 m 0. s 80, Deck = -52°26'25."1 
using the XMMSAS tool barycen. Finally, we selected 
counts within the energy range 0.4—2.5 keV to maxi- 
mize the signal-to-noise ratio of the pulsed signal before 
beginning our timing analysis. The numbers of selected 
counts are given in Table 1 (background was estimated 
to contribute less than 15% in each dataset). 



3 Phase-Coherent Timing Analysis 

Phase-coherent timing analysis requires careful spacing 
of individual observations such that an extrapolation of 
the measured phase model, <f)(t) = (f>(t n ) + v(t — t ) + 
\v(t — t a ) 2 + for a given observation or set of ob- 
servations is precise enough to predict the phase to the 
next observation to much better than a pulse cycle. The 
advantage of this approach is that one can achieve far 
more precise measurements of the pulse frequency and 
higher derivatives than by using independent pulse fre- 
quency measurements with the same total exposure. This 
approach is commonly applied to all types of pulsars, in- 
cluding Anomalous X-ray Pulsars (e.g. Gavriil & Kaspi 
2004), Soft Gamma Repeaters (e.g. Woods et al. 2002), 
and radio-quiet Isolated Neutron Stars (Kaplan & van 
Kerkwijk 2005a,b). 

3.1 Pulse Phase Fitting Technique 

As the frequency error (Sis) in an individual observation 
is inversely proportional to its duration, 8vj oc TJ , the 
longer central exposure served as our reference point. 
We measured the pulse frequency during this observa- 
tion first via a Z 2 search (see §3.2) and then refined this 
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Table 1 XMM-Newton observation log for IE 1207.4-5209. 



# 


ObsID 


Date 


Central epoch 


Span 


PN Exp. a 


r b 


Counts c 


max Z\ 








(MJD) 


(ks) 


(ks) 


(arcsec) 






1 


0304531501 


2005 Jun 22 


53543.600542 


15.1 


10.6 


35 


13,233 


41.1 


2 


0304531601 


2005 Jul 05 


53556.141927 


18.2 


12.7 


35 


12,858 


37.0 


3 


0304531701 


2005 Jul 10 


53561.404084 


20.5 


14.3 


20 


17,651 


43.1 


4 


0304531801 


2005 Jul 11 


53562.456311 


63.4 


44.4 


35 


56,804 


120.4 


5 


0304531901 


2005 Jul 12 


53563.335586 


9.6 


6.7 


20 


8,559 


19.1 


6 


0304532001 


2005 Jul 17 


53568.112822 


16.5 


11.5 


35 


14,696 


81.2 


7 


0304532101 


2005 Jul 31 


53582.691485 


17.7 


12.4 


20 


15,524 


26.5 




Sum 






161.0 


112.6 




139,325 





a Effective source exposure times after filtering. See text for details. 
b Source extraction radius used for event selection. 
c Number of counts used for timing analysis. 



measurement as follows. We split the observation into 4 
segments and folded these segments on the measured fre- 
quency to generate pulse profiles for each segment. Next, 
we cross-correlated each pulse profile with a high signal- 
to-noise pulse template and measured phase offsets. The 
pulse template is first derived from the central observa- 
tion folded at the initial frequency. The phase offsets for 
the 4 segments were fitted to a straight line and the slope 
of this line was added to the initial frequency to produce 
our refined frequency. The short gaps between the cen- 
tral exposure and the adjacent exposures were expected 
to preserve the phase information, i.e. the propagated 
phase error (e.g., between the 4-th and 5-th observations, 
o~<fi = 5v (t$ — ti)) was expected to be <C 1 cycle, which 
would mean that no pulse cycles are missed in the phase 
model. As one incorporates more and more data over 
a wider time span, the precision of the phase model im- 
proves, and one can tolerate larger gaps between observa- 
tions. Note that the template pulse profile is updated as 
more data are included until the full data set is utilized. 
By the time we incorporated the measured phases from 
the first and final observations into our fit, it became 
clear that the phase offsets did not conform to a simple 
linear trend, and a quadratic term (oc v) was added to 
the phase model, <f>{t) = <fi(to) + v(t — to) + \v{t — to) 2 . 
However, even the inclusion of the quadratic term did 
not reduce the variance of the phase residuals to the 
point where we obtained an acceptable fit (x 2 = 19.2 for 
8 degrees of freedom; see Fig. 1). 

The poor fit to the quadratic phase model indicated 
that we either converged on an alias solution or 1E1207 
exhibits significant "phase noise" 1 on a time scale of 
weeks. An alias timing solution would be when there 
are an incorrect number of cycle counts between consec- 
utive observations. Phase noise can be characterized in 
many ways such as the presence of a strong cubic term 

1 In this context, we refer to phase noise simply as devia- 
tions from our simple quadratic phase model beyond statisti- 
cal errors. Note that an alternative definition of phase noise 
has specific meaning in the context of pulsar timing noise 
(e.g. Cordes & Helfand 1980). 




-20 -10 10 20 

Time (Days since MJD 53563.148) 

Fig. 1 Pulse phase residuals from XMM-Newton observa- 
tions of 1E1207 during the 2005 observing campaign for 
model MODI (Table 2). Top: Phase residuals minus a lin- 
ear trend. Bottom: Phase residuals minus a quadratic trend. 
Note that the central (longest) observation is split into four 
segments of equal spans. 

(oc i>), white noise, periodic variations, etc. To ensure 
that the poor \ 2 value in the solution we found is not the 
consequence of misidentified cycle counts, we employed 
a technique used for timing noisy rotators such as Soft 
Gamma Repeaters (Woods et al. 2006). In this technique, 
we measure the pulse phase and frequency at each of the 
7 observing epochs. The phase for each observation was 
measured by folding the data from each observation on a 
pulse ephemeris of constant frequency determined by the 
central observation and computing the phase difference 
between this profile and a template pulse profile. The 
pulse frequencies for the short observations were mea- 
sured by splitting the data into two segments of equal 
duration, folding each segment on the pulse frequency 
measured for the central observation, measuring phase 
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shifts for each pulse profile relative to the pulse tem- 
plate, and fitting these two pulse phases to a line to de- 
termine the local pulse frequency. Finally, we perform a 
least-squares fit to the set of 7 phases and 7 frequencies, 
where we vary the number of cycles between consecutive 
observations by integer increments (Woods et al. 2006). 
This provides a family of solutions to the full data set 
which define the pulse phase evolution according to a 
quadratic model covering the 40-day time interval. None 
of the solutions (for a quadratic phase model) provide 
a statistically acceptable fit to the data. For all possi- 
ble solutions, the "null hypothesis probability" (i.e. the 
probability of measuring the large \ 2 values by chance, 
assuming that the model is correct) is very small. Fit pa- 
rameters for the top three solutions (ranked in order of 
increasing x 2 ) are given in Tabled We chose to consider 
a limited number of solutions; therefore, we selected only 
the solutions that had a probability of getting the mea- 
sured x 2 by chance of 10~ 5 or larger (three solutions). We 
found that the best-fit model is equivalent to the solution 
we identified via our bootstrap phase-fitting method de- 
scribed earlier. Phase residuals for this model are shown 
in Figure^ Since none of the identified solutions provide 
a statistically-acceptable fit to the data, we conclude that 
lEf 207 does, in fact, exhibit significant phase noise on 
a time scale of weeks. It seems unlikely that the excess 
noise we observed is due to underestimating our phase er- 
rors. Analysis of XMM-Newton data from other pulsars 
using the same software, although covering time spans 
shorter than 40 days, have consistently yielded reduced 
X 2 values of ~1 (e.g. Woods et al. 2004). 

The presence of the phase noise does not allow us to 
unambiguously phase-connect the complete data set and 
thus measure a unique frequency and frequency deriva- 
tive for the full 40-day time span. To place some con- 
straints on the frequency derivative during our observing 
sequence, we employed a Monte-Carlo simulation to esti- 
mate statistical significance of the multiple solutions. For 
this simulation, we first had to choose a model for the 
phase noise. We selected two models: (i) a cubic phase 
term and (ii) white noise. In both cases, the amplitude 
of the model noise variance was equal to the total vari- 
ance in the top three fits to the data minus the statis- 
tical variance. In our simulation, we generated phases 
for each observing epoch which included three compo- 
nents: the model phase (including v and v terms), Gaus- 
sian measurement noise, and the model phase noise. In 
addition, we simulated frequency measurements at each 
epoch assuming Gaussian measurement noise (i.e. we ne- 
glect phase noise on the time scale of the observation 
duration). For each phase model, we generated 10 5 real- 
izations and fit for the cycle counts between consecutive 
epochs as we did for the measured data to identify all 
possible timing solutions for each realization. In each re- 
alization, we identified the rank of the true timing solu- 
tion in terms of x 2 . The most constraining results were 
obtained from the white noise model for the phase noise. 



For this model, we found that the true timing solution 
was among the top three solutions (ranked in order of 
X 2 ) 90% of the time and was the top solution 65% of the 
time. Assuming white phase noise, our simulation sug- 
gests that we can be 90% confident that the true pulse 
ephemeris for 1E1207 is MODI, MOD2 or MOD3 given 
in Table Similarly, these results suggest there is a 65% 
chance that MODI defines the appropriate cycle counts 
between observing epochs, and hence, reflects the correct 
pulse ephemeris. 

The differences between the three pulse ephemerides 
listed in Table amount to small differences in the cy- 
cle counts between the four outer observations in our 
observing sequence (i.e. a few additional or less cycles 
between observations 1 and 2, 2 and 3, 5 and 6, and 6 
and 7). In fact, we can only be sure of the cycle count 
accuracy between the three central observations (obser- 
vations 3, 4 and 5 in Table 0). To show this explicitly, 
we fit for the cycle counts between the central three ob- 
servations as we did for the full data set, only we limited 
the order of the phase model to be first order on ac- 
count of the short time span (2 days). We measure a 
difference in x 2 of 91 for 3 degrees of freedom between 
the best-fit ephemeris identified in our search and the 
next closest. Clearly, we were able to phase-connect this 
subset of the data and unambiguously identify the local 
pulse frequency (v — 2.35776187(31) Hz over the time 
range 53561.328 to 53563.347 MJD TDB). 

Although the method described in this section is very 
efficient, it has some limitations. For large cycle count 
corrections between consecutive observations, the local 
pulse ephemeris will change considerably as will the folded 
pulse profile. In turn, the pulse phase measurement will 
likely also be affected. In practice, the differences in the 
pulse shapes of 1E1207 for the three pulse ephemerides 
reported here are insignificant. For very large cycle count 
corrections, where this effect becomes important, the x 2 
contribution from the frequency measurements begin to 
dominate the total x 2 , and these peaks are effectively 
suppressed. Even so, this method is relatively new and 
not extensively tested. To verify the results obtained 
with this technique, we employ the Z 2 test, a traditional 
approach to X-ray timing. 



3.2 The Z 2 test 

The Z 2 statistic (e.g., Buccheri et al. 1983) is defined as 
follows: 



Z 2 = 



N 



E 

fe=i 



N 



N 



cos 2-nk(j)i ) + I sin 2-Kkcj>i 



-(1) 



where 4>i — v{ti — to) + v(U — to) 2 /2 + ... is the phase of 
i-th event, ti — to is the event arrival time counted from 
an epoch to of zero phase, n is the number of harmonics 
involved in the test, and N is the number of events. For 
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Table 2 Candidate pulse ephemerides for IE 1207.4-5209 for 2005 June through July. 



Model 


Epoch 


u b 


i> 


X 2 /dof 


Null Hypothesis 


zl 


zl 




(MJD TDB) 


(Hz) 


(10~ 12 Hz s- 1 ) 




Probability 






MODI 


53563.148 


2.357761722(16) 


+0.890(22) 


25.2/11 


2.9 x 10~ a 


323.0 


333.5 


MOD2 


53563.148 


2.357762311(16) 


-2.651(23) 


32.6/11 


2.3 x 10 -4 


314.9 


325.8 


MOD3 


53563.148 


2.357761720(17) 


+1.607(24) 


37.6/11 


3.6 x 10" 5 


318.5 


325.2 



a Pulse ephemerides are valid over the time range 53543.547 to 53582.746 MJD TDB. 

6 Numbers given in parentheses indicate the la error in the least significant digit(s). The statistical errors are inflated by a 
factor (x 2 /dof) 1/2 . 



a signal with a nearly sinusoidal pulse profile, such as 
observed from 1E1207, the Z\ (Rayleigh) test is known 
to give excellent results. For a sinusoidal signal, the ex- 
pected peak value of Z\ is -/V/ 2 /2, where f p is the pulsed 
fraction. For 1E1207, the pulsed fraction was measured 
to be f p =8%-12% (Zavlin et al. 2000; Pavlov et al. 
2002b). The peak Z\ values found in the individual data 
sets (Table 1) are in a reasonable agreement with those 
given by this estimate. 

The Z^ test has been used for a phase-coherent tim- 
ing analysis of several observations spread over a large 
time span by Mattox et al. (1996) and Zavlin et al. 
(1999), and we follow the approach described by those 
authors. To account for the phase connection, we apply 
the Z^ test (for n = 1 and 2) to the whole data set of 
seven observations. To determine the parameters v and 
v of the quadratic phase model, we calculated the Z 2 
on a dense two-dimensional grid [v — 2.3577 Hz = 41- 
81 /iHz, \v\ < 1 x 10~ n Hz s _1 ], with v and i> spac- 
ings of 0.02 /xHz and 2 x 10~ 14 Hz s _1 , respectively. A 
contour map obtained with the Z\ statistic is shown in 
Figure [21 Because of the cycle-count ambiguities dur- 
ing the gaps between the consecutive observations, the 
map shows multiple peaks, one of them corresponding to 
the true v^v solution and the others being aliases. The 
first, third and fourth highest peaks in this map corre- 
spond to MODI, MOD3 and MOD2, respectively (see 
Table 2). The top three peaks in a similar Z\ map are 
at the same v,v as MODI, MOD2 and MOD3, respec- 
tively. If the phase connection between separate data sets 
were perfect, then the peak corresponding to the true 
solution would be much higher than the aliases. How- 
ever, in our case the difference between the heights of 
the peaks turned out to be too small to single out a 
unique solution. For instance, in addition to the highest 
peak in the Z\ map, Z\ max = 323.0 at v = 2, 357, 761.72 
/iHz, v — +0.90 x 10 -12 Hz s _1 , we see four peaks with 
310 < Z\ < 320 in FigureEl at different v,v values. Sim- 
ilar to the method described in §3.1, the differences in 
peak values of v,i> correspond to different (integer) num- 
bers of cycles (~8x 10 6 ) during the full observational 
time span T = 3393.8 ks. We are not aware of statistical 
criteria to estimate significance of separate peaks in this 
approach, and we can only assume that the solutions cor- 
responding to several highest Z\ peaks cannot be ruled 
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17-2.3577 [^Hz] 

Fig. 2 Contour plot of the Z\ power on the v-v grid. The 
purple, green and red contours correspond to Z\ = 250, 280, 
and 310, respectively. 



out on statistical grounds. We also note that the lack 
of perfect phase-coherence is supported by the fact that 
the largest Z\ is much smaller than Ylj=i — 346.8 
(at the same v, i>), the value we would expect to obtain 
for perfect phase connection. Thus, the results of the Z^ 
search for the 1E1207 frequency and frequency deriva- 
tive are generally consistent with the results reported in 
63.1. 



4 Origin of the erratic spin behavior 

The deviations from monotonic spin-down in 1E1207 are 
substantial, and they manifest on timescales of years to 
as short as possibly weeks as evidenced by the phase 
noise detected here. We now consider four possibilities 
for both the erratic long-term spin behavior and short- 
term phase noise in 1E1207: (i) frequent glitching, (u) 
accretion and propeller torques from a circumstellar de- 
bris disk, (Hi) timing noise in an isolated neutron star, 
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and (iv) orbital Dopplcr shifts caused by the presence of 
a binary companion. 

For any glitch model, the glitch frequency and am- 
plitude would have to be very high to account for the 
observed spin variability (see Zavlin et al. 2004). More- 
over, the most viable timing solutions indicating spin 
down over the 40-day observing span differ from the 
long-term spin-down of 1E1207 over the last 5.5 years 
(~ —4 x 10~ 14 Hz s _1 ) by an order of magnitude. For 
example, the most likely spin-down solution (MOD2) has 
a frequency derivative more than one order of magni- 
tude larger. Because only a very contrived glitch model 
could account for the long-term frequency history, and 
this model would provide no explanation for the short- 
term phase noise, the glitch model is effectively excluded 
by these observations. 

Debris disks left over from the supernova explosions 
that produce neutron stars could alter the spin evolu- 
tion of the central neutron star via accretion and pro- 
peller torques (Zavlin et al. 2004). If the spin-up rate 
of 1E1207 during the 40-day interval were equal to the 
values measured for MODI or MOD3, then the mass ac- 
cretion rate would have to be very large (to > 3 x 10 16 
g s _1 ). Such a large accretion rate would require a large 
increase in X-ray luminosity which is not observed. Even 
the spin-down solutions would require significant optical 
and IR emission from the disk. Deep IR observations of 
1E1207 have shown no indication of even a cool, passive 
debris disk (Wang, Kaplan & Chakrabarty 2006). Thus, 
it appears unlikely that a debris disk is the cause of the 
spin variability in 1E1207. 

Timing noise (irregular evolution of the pulse phase 
with time) is a ubiquitous phenomenon in isolated neu- 
tron stars. This variability is in addition to the usual vari- 
ation caused by magnetic braking. It has been demon- 
strated that the magnitude of these irregular variations 
depends upon the spin-down rate of the pulsar (e.g., 
Cordes & Helfand 1980). Millisecond pulsars show the 
smallest timing noise while magnetars exhibit very strong 
timing noise. In the case of magnetars, these variations 
manifest as changes in the effective spin-down rate of 
up to factors of 5 on a time scale of years. For a conve- 
nient (albeit crude) description of timing noise, Arzou- 
manian et al. (1994) introduced a "stability parameter" 
defined by the following equation: A\ ogt — log(|i>|i 3 /6^), 
where t is the time during which the pulse phase has 
been monitored (t — 10 s s is a commonly used charac- 
teristic time), and v is the formal value of the second 
frequency derivative obtained from fitting a cubic model 
to pulse phases (it is much larger in magnitude than the 
actual v for noisy pulsars). Third-order polynomial fits 
to the 1E1207 phase residuals of the top three candidate 
timing solutions yielded insignificant measurements of v. 
The timing observations of 1E1207 during 5.5 years were 
too sparse to fit the pulse phases with any model. There- 
fore, to estimate v and A%, we fitted the dependence 
v(t) = vq + i>(t — to) + i)(t — to) 2 /2 to the frequency his- 
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Fig. 3 The timing noise parameter As after Arzoumanian 
et al. (1994) for 119 radio pulsars, the isolated neutron 
star RX J0720.4-3125 (Kaplan & van Kerkwijk 2005b), 
four Anomalous X-ray Pulsars (blue squares - Gavriil & 
Kaspi 2004; blue plus signs - Heyl & Hernquist 1999), two 
Soft Gamma Repeaters (Woods et al. 2002) and the CCO 
IE 1207.4-5209. 



tory covering the last 5.5 years. (The same exercise was 
performed by Heyl & Hernquist 1999 for some Anoma- 
lous X-ray Pulsars, and it was shown to be reasonably 
accurate by Gavriil & Kaspi 2004.) Choosing to — 52700 
MJD TDB, we found v Q = 2, 357, 762.8 ± 0.2 /zHz, v = 
(-3.4±0.6)xl0~ 14 Hz standi/ = (5.9±1.9) x 10~ 22 Hz 
s~ 2 , which translates to a timing noise level of As = 1.6. 
In FigureOl we show the period derivative versus the tim- 
ing noise parameter A$ for 126 isolated pulsars of various 
flavors as well as for the CCO 1E1207. The isolated pul- 
sars fall along a relatively well-defined locus, from the 
quiet millisecond pulsars to the noisy Soft Gamma Re- 
peaters, while 1E1207 stands out from this trend with 
an anomalously large timing noise strength, some 2—4 
orders of magnitude higher than isolated pulsars at sim- 
ilar spin-down rates. Although such an estimate for the 
timing noise parameter is, by necessity, very crude, its 
enormously high magnitude, together with the gross in- 
consistency of the local (June-July 2005) spin-down rate 
with the long-term average, suggest that the erratic fre- 
quency behavior in this source is not due to the same 
effect that causes timing noise in other isolated neutron 
stars. 

The most straightforward explanation for the long- 
term spin variations in 1E1207 is the presence of a bi- 
nary companion. Note that this model cannot explain 
the observed short-term phase noise. Current IR and 
optical limits for 1E1207 exclude main sequence com- 
panions earlier than M5 and even white dwarfs with ef- 
fective temperatures greater than ~10 4 K (Fesen et al. 
2006; Wang et al. 2006). Allowable companion masses 
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Fig. 4 Candidate orbital solutions for 1E1207 consistent 
with our MOD2 (top) timing solution and MODI (bottom) 
timing solution. 



5 Conclusions 

We observed 1E1207 with XMM-Newton seven times 
during the course of a 40 day interval in an effort to mea- 
sure the local pulse frequency and frequency derivative 
with high precision. Due to unanticipated phase noise, we 
were unable to phase-connect the full data set. From sys- 
tematic pulse ephemeris searches, we identified a number 
of possible timing solutions, none of which had spin-down 
rates close to the average spin-down rate of 1E1207 over 
the last 5 years. The interpretation of the erratic long- 
term timing behavior of 1E1207 in terms of the usual pul- 
sar timing noise would require timing noise levels much 
higher than seen in isolated neutron stars with compara- 
ble spin-down rates. The most plausible explanation for 
the erratic long-term pulse frequency evolution of 1E1207 
is the presence of a binary companion, although the cur- 
rent data do not allow us to place strong constraints 
on system parameters. Further efforts at phase-coherent 
pulse timing observations of 1E1207 are required to (i) 
unambiguously identify the nature of the long-term pulse 
frequency variations, and (ii) confirm and further inves- 
tigate the observed short-term phase noise. Just one ad- 
ditional 40-day observing sequence with higher sampling 
would yield far better constraints on orbital parameters 
should 1E1207, in fact, possess a binary companion. 
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are less than 0.2 M Q for late-type stars (Moody et al. 
2006, in preparation). Even with such low-mass com- 
panions, the resulting Doppler shifts are large enough 
to account for the frequency variations in 1E1207. Us- 
ing the archival spin frequencies in combination with 
the phases from our three candidate timing solutions, 
we fit the data to a circular orbital model whose phase 
evolution is defined by the following equation: cf){t) = 
4>(t ) + u(t - t ) + \v\t - t ) 2 +Asmu(t - £ ). This is 
the same equation as given in §3.1 with an additional 
sinusoidal term to account for the orbital Doppler shifts. 
We identified a family of acceptable orbits for each of the 
three timing solutions listed in Table HJ The full set of 
allowable timing solutions are too numerous to list. We 
can place only very crude constraints on the orbital peri- 
ods to fall between 120 and 600 days. The mass functions 
range between 1 x 10~ 7 and 5 x 10 -5 M for acceptable 
orbital solutions. For a 90° inclination and a 1.4 M 
neutron star, the corresponding companion mass range 
is 0.007 to 0.05 M Q , well within the existing limits on 
companion masses from IR and optical observations. For 
illustrative purposes, we show two example orbital solu- 
tions that are consistent with the existing timing data 
for 1E1207 (Figure EJ. 
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